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A  I'lACHINE  PROGRAM  FOR 

COMPUTING  NONLINEAR  VISCOELASTIC  WAVE  PROPAGATION  IN  SOILS 

ABSTRACT 

A  numerical  method  and  its  machine  program  for  calculating  nonlinear  viscoe¬ 
lastic  wave  propagation  in  soils  are  described.  Solutions  give  particle  velocity; 
stress  and  strain  as  functions  of  time  and  distance  from  the  stress  source. 
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I.  INTRODUCTION 


This  report  is  primarily  a  reference  to  the  BRLESC  (Ballistic  Research  Lah- 
atories’  Electronic  Scientific  Computer)  computational  procedure  used  at  the  Bal¬ 
listic  Research  Laboratories  to  solve  a  system  of  nonlinear  differential  equations 
describing  the  viscoelastic  propagation  of  stress  waves  in  soils. 

Studies  of  this  nature  have  been  pursued  in  the  past.  Seismic  forces  and 

their  associated  motion  have  been  measured  under  known  conditions  to  check  theo- 

12  5* 

retical  development  ’  ’  .  Simulation  techniques  have  been  used  for  the  longer 

duration  pulses  produced  by  nuclear  explosions.  Measurements  in  the  free  field 
have  shown  that  more  consideration  should  be  given  to  the  phenomena  occurring 
after  the  first  stress  peak.  In  order  to  evaluate  theory  on  the  diffraction  of 
stress  waves  around  an  enclosirre,  a  simple  structure  has  been  buried  in  a  soil 
filled  shock  tube.  The  techniques  are  being  improved  for  measurements  of  the 
required  parameters  associated  with  stress  propagation. 

When  a  nuclear  device  is  detonated  in  air  near  the  ground,  a  large  area  of 
the  ground  under  the  point  of  detonation  is  loaded  nearly  simultaneously.  Under 
such  loading,  a  uniform  or  horizontally  layered  half  space  woiild  experience  motion 
nearly  perpendicular  to  the  surface.  Plane  wave  theory  can  be  used  to  predict 
with  sufficient  accuracy  the  ground  parameters  near  the  surface  which  is  close  to 
the  point  under  the  detonation. 

Since  in  the  high  stress  region,  we  cannot  use  available  linear  elastic  theory, 
nonlinear  stress-strain  assumptions  must  be  made.  Mercado  assumed  a  nonlinear 

viscoelastic  theory  (time  dependent)  which  has  been  shown  to  be  present  under  some 
extreme  conditions.  The  quantitative  importance  of  time  dependency  has  not  yet 
been  established  for  the  phenomena  of  interest  here,  and  this  effect  may  not  be  of 
major  concern  for  many  materials. 

This  work  is  a  theoretical  phase  of  the  contract  on  propagation  of  under¬ 
ground  shock  waves  contracted  to  Rensselaer  Polytechnic  Institute  by  the  Terminal 
Ballistic  Laboratory  at  Ballistic  Research  Laboratories. 

Information  obtained  from  such  stress-wave  studies  are  required  for  free  field 
predictions  and  will  be  useful  in  the  development  of  techniques  for  destroying  or 
protecting  underground  structiires  from  nuclear  blast. 


* 

Superscript  numbers  denote  references  listed  at  the  end  of  this  report. 
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II.  THEORY 


4^6 

The  theory  discussed  here  is  that  of  Mercado  ’  ^  which  is  based  on  Malvern's 
V  3 

theory  for  stress  propagation  in  strain-rate  sensitive  materials. 


An  analytical  representation  of  the  stress- strain  curve  is  obtained  by  fitting 

Q 

an  empirical  formula  suggested  by  Osgood  of  the  form  (Reference  U,  Equation  2) 
(Il-l)  €  =  I  +  R(|)^ 

to  the  static  data.  In  the  above  equation  e  =  strain^  a  =  stress,  E  =  Young's 
modulus  and  R  and  S  are  static  constants  of  a  material.  With  this  law,  the  basic 
constitutive  eqiiatlon  defining  the  material  (Reference  4,  Equation  lO)  takes  the 
form 


(II-2) 


b 

a  ^t 


1 

a 


y 


where  t  represents  time,  a  and  b  are  coefficients  in  the  viscoelastic  law,  and  we 
choose 

(II-3)  e(c)  =  I  +  K(|)  , 

where  H, K  are  parameters  in  the  stress-strain  law.  Mercado  states  that  there  is 
no  experimental  data  known  from  which  the  dynamic  constants  a  and  b  may  be  directly 
determined^^ .  It  is  known  from  theory  that  b/a  >  E,  therefore,  b/a  is  arbitrarily 
chosen.  The  lowest  order  approximation  to  the  dynamic  stress-strain  law  is  obtained 
by  truncating  the  time-wise  expansions  of  a  and  e  after  the  first  derivative,  giving 


(II-4)  6(cr)  +  a  =  f(e)  +  be^, 

where  we  arbitrarily  choose 


(II-5) 


f(€)  =  e 


1  +  ge 
1  +  ae^ 


a  and  p  being  parameters  in  the  stress-strain  law. 
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The  dynamic  "behavior  of  materials  is  obtained  by  simultaneous  solution  of  the 
dynamic  stress- strain  law  and  the  laws  for  the  conservation  of  mass  and  m^omentum. 


lor  a  first  approximation  to  the  dynamic  behavior  of  soils,  we  use  the  stress- 
strain  lav  the  equation  for  the  conservation  of  mass 


(II-6) 


e.  =  V  , 
t  X-’ 


and  the  equation  for  conservation  of  momentum 


(II-T) 


Here  x  is  Lagrangian  distance,  v  is  material  velocity,  and  p  is  density.  Linear 
combinations  of  Equations  II-4,  II-6  and  II-7  give  the  following  set  of  charac¬ 
teristic  equations.  Along  ^  =  0  (particle  path), 

CL  L 


(II-8) 


do  -  —  de  =  -—  g(a)  -  f(e)  dt; 
a  a  1_ 


=  -F(cr,e)dt 


and  along  ^ 
where 


(positive  cnaracterlstic). 


c  =  /£H 
■v'  a 


F  =  i  g(a) 

cl 


(II-9) 


do  -  c  dv  =  -—  g(cr)  -  f(e)  dt 

^  L 


=  -F(cr,£)dt; 


dx 

and  along  ^  = 


(II-IO) 


-  (negative  characteristic), 

Q.C 


+  C  dv  =  -i  g(a) 


f(G)  dt 


=  -F(a, G)dt. 


These  equations  are  supplemented  by  the  jump  conditions  for  a  discontinuous 
wave  propagating  into  a  stationary  medi-um  in  which  a  =  0,  e  =  0  and  v  =  0.  By 
replacing  e,  and  v  respectively  in  (II-6),  (lI-7)  by  aH(|),  and 

vH( g ),  where  H(  § )  =0  for  |  <  0  and  H(  O  =  1  for  I  >  0  and  =  Dirac  6,  we 

integrate  the  results  with  respect  to  |  from  -i  to  and  letting  |  0  we  obtain 

the  jump  relationship 


f 

V. 


II-ll) 


a  = 


b 

-CV  =  —  €. 
a 


These  equations  represent  the  conservation  of  mass  and  momentum  across  a  disconti- 
uity  moving  into  an  landisturbed  medium  with  the  velocity 


(11-12) 


dx  b 
dt  ac' 


which  is  also  a  positive  characteristic. 

III.  NUMERICAL  FORMULATION 

To  compute  the  stress,  strain,  and  particle  velocity  behind  the  wave  front, 
we  integrate  the  conservation  equations  using  the  known  values  along  the  wave  front 
and  the  initial  and  bo\nidary  values. 

The  characteristic  curves  are  used  as  the  coordinate  system  to  extend  the 
solution  into  the  domain  of  disturbance. 

Using  a  grid  of  characteristics  in  the  t,x  plane  (Figure  l)  a,  e,  and  v  are 
evaluated  at  each  point  (m,n). 

Along  the  wave  front  (m  =  O),  (lI-9)  and  (ll-ll)  hold.  Differentiating  Equa¬ 
tion  II-ll  with  respect  to  t  and  substituting  into  Equation  II-9>  we  obtain 


(III-I) 


H  =  |^g(a)  -  f(€)  =  -|F(0,e). 
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FIGURE  1. 
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In  difference  form  this  eq.uation  becomes  (see  Figure  2) 


(III-2) 

or 

(III-3) 

With  this 

(III-1+) 

(III- 5) 
Initially 

and 


FIGURE  2. 


-  cr 

o,n+l  o,n 
2 


(1-1)  (i-l)\ 

.%,n+l^  ^o^n+ly  ^  ^^‘^Ojn^  ^o,n^ 


(i) 


~  A  ^  ,,]  +  F(cr  ,e  ) 

o,n+l  o,n  8  \  o,a+l'  o,n+l/  '  o,n  o,n 


/(i-1)  (1-1)\ 


value  of  a  at  (o,n+l),  we  have  e  and  v  from  (ll-ll), 


(i)  ^  (i) 

^o,n+l  ~  b  ^o,n+l' 


(i)  1  (1) 

V  =  — ~  O'  • 

o,n+l  c  o,n+l 


a  =  a  ,  e  =  ^  o 
0,0  o'  0,0  b  o 


0,0 


’c"  ^o ' 
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The  set  (III-3),  (III-5)  are  iterated  over  i  at  each  point  (o,n), 

until^  for  some  arbitrarily  small  positive  number  5, 


(III-6) 


rro) — 


<  6? 


n  is  then  increased  to  some  maximum  value.  Thus,  the  m  =  0  characteristic  is 
obtained. 


For  the  m  =  1  characteristic,  the  first  point  to  be  computed  is  the  boundary 

point  on  the  line  x  =  0.  On  this  line  we  assume  some  constant  value  a  for  the 

o 

stress.  Then  the  general  boundary  point  (m,m)  (see  Figure  3),  we  compute  as  fol¬ 
lows  : 


The  time  and  Lagrange  distance  are  obtained  from 


(III-7) 


t 

m,m 


ir-l,m-l 


+  At, 


(III-8) 


X 

ra,ra 


0. 


The  stress  is  the  specified  value 


(III-9) 

II 

B 

a  . 
o 

Let 

(III-IO) 

m-l,m-l  m-l,m 

m,m 

2  " 

(Ill-ll) 

v^°^  = 

V  +  V 

m-l,m-l  m-l,m 

m,m 

2 
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The  difference  form  of  (lI-8)  gives  for  the  strain 


(III-12) 


.(i)  _ 


m,m 


f{(j  ,  e  T  )  +  F(  ( 
^  “h - - 2 - 


X--L^ 


a  ,  e 
0^  m,m 


and  the  difference  form  of  (lI-9)  in  which  da  =  0  gives  for  the  velocity 


(III-13) 


m,m 


1  / 

V  , - (a 

m-l,m  c  o 


At 

2c 


This  process  is  repeated  until  for  some  i,  the  convergence  test 


(III-14) 


•  rn  m  *  y 


m.m 


0. 


ITT) 

m,m 


_num 


5 


is  satisfied,  where  0  =  €  and  v,  and  where  if  the  denominator  is  zero,  this  con¬ 
vergence  test  is  bypassed. 


With  this  boundary  point  computed,  we  progress  up  along 
acteristic.  Each  point  along  this  characteristic,  say  point 


the  positive  char- 
(m,n)  (see  Figure  i+). 


t 


FIGURE  4. 
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Tlie  time  and  Lagrange  distance  are  obtained  from 


(Ill-lf) 


=  t 


m,n-l 


+ 


At 

2^ 


(III-16) 


X 

m,  n 


b  /Atv 

^in,n-l  ac  2'’ 


L-t 


(III-IT) 


(o)  _  ^m-l,n 


m^n-l 


m,  n 


(iii-ie) 


(o)  Si,n-1 

^m,n  "  2 


(III- 19) 


f  \  v  ,  +  V  T 

(o)  Tn-l,n  m,n-l 

\,n  ■  2 


The  difference  forms  of  (lI-8),  (lI-9),  and  (ll-io)  solved  simultaneously  give 


(III-20) 


(i)  m,n-l  m-l,n  c, 

a  =  — ^ - o - ^ - o(v  ^ 

m^n  2  2  ra^n-l 


At  F^^m,n-1^ 

■^L 


e 

m,  n- 


1 


)  +  F 

~2 


m-l,n^  m-l^n 


' _ \  m,  n^  m,n 
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(III-21) 


a 


m,n  m-l,  n-1 


(i)  ^ 

O  -  CT  ,  ^  J 

,  m,n  m-l,n-iy 


(III-22) 


F(0  1^1  i)  + 

At  m-l^n-l,  m-l,n-l'^  \  m,n^  m,n  / 

B  2  " 


(±) 

f  •  \  o'"  '  -  a  ^ 

=  V  + 

m,  n  m,  n-1  c 


At 


F(0  ,  e  -I ) 

m,n-l  rti,n-I 


^{i-l)^(l-l) 
m,n-^  mjn 


This  process  is  repeated  for  each  i  until  the  convergence  criterion 


(m-23) 


^m,n 


0, 


(1-iL 


A 


m,n 


m^n 


< 


5 


is  satisfied  for  0=0,  e,  and  v. 


IV.  ORGANIZATION  OF  CALCULATION  AND  PROGRAM 

For  each  m  characteristic  that  we  compute,  we  require  values  on  the  m-1 
characteristic;  therefore,  the  memory  size  limits  the  number  of  points  that  can  be 
computed  on  each  characteristic.  In  order  not  to  limit  the  domain  of  computation, 
it  was  divided  into  zones.  For  this  purpose,  four  codes  are  written,  differing 
only  in  the  input  arrangements.  Tnese  foior  codes  can  be  used  to  cover  the  entire 
domain  (Figure  5)- 
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All  option  provider,  a  choice  on  the  nunber  of  computed  points  to  be  printed. 
On  the  boundaries  of  each  code  zone,  however,  values  at  all  points  are  printed, 
.'.inec  they  may  be  the  input;  data  for  sulueciuent  zones.  Code  1  would  suffice  to 
.;o;;irlcte  the  computation  if  the  memory  were-  sufficiently  large  to  accommodate  all 
the  -points  on  an;-  given  characterist;ic . 

Tlie  compiler  language  used  is  the 

CODE  SllIBOLS 


CODE  NAME 

DEFINITIONS 

A 

1 

AIjPHA 

a 

a 

AN3 

Printing  freq'iency 

B 

(1,2, 3, etc.  ) 
b 

BEE 

a 

1- 

( 

BEGIN 

r+  no  shock  at  wave  front 

!  EFl'A 

L-  shock  at  wave  front 

P 

c 

c 

DELE 

At 

DELTA 

5 

1  ^ 

e  (strain) 

i  EG 

E  (Young's  modulus) 

i  H 

H 

K 

K 

M 

m 

N 

n 

RHO 

P 

SIGMA 

a 

SIGO 

^0  TEOHNIOAL  LIERASI 

T 

t  BLP  cl3 

V 

abeedeen  proving  ground,  m 

■V  STEAP-TL 

X 

X 
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CODE  1 

NUMERICAL  RESULTS  FOR  THEORETICAL 
VISCOELASTIC  WAVE  PROPAGATION  FOR 
SHOCK  STUDIES. 

NON-LINEAR 

GROUND 

/startX 

I  2i 

Gi 

START 

print  ^  N 

IS  LESS  THAN  M  > 

ywl  i' 

(  VPI-  VP7/3  ) 

(  8I0PI  -8I0P7/5) 

PI  FORM  (9-10)  1-7)  5-IMO) 
LAST  (EPI)  VPI)  SIOPI) 


2  >1 

SI(}PI,  F>  SlOO 
EPI.  F*  I/e  X  SICO 
VPI,  F>l/C  X  SlOO 
PRINT-  FORMAT  (FI)-(MC)NC) 
T)X)EPI.F)VPI,F)SIQPI  F) 
SET  (3.7ICX*  634) 


INC  (Ft  F  «-S  ) 


SlOPliF-  3)*  0^ 


OsSlG  PI,  (F-3  )/E0  ♦ 
K•EKP(H•LOO 
(SIOPI.(F-3)/EO)) 


PRINT  labeling 
PRINT  BLANKS 


Z«ANS 
B>  BEE«A 
SET (F*0) 

C=  SORT(M) 

MC*0 
NC  s  0 
T-.  0 
X  *  0 

PRI  NT  -  FOR  MAT  ( F I )  •  ( SlOO) 
BEGIN)ANS)B)  C) 
LABEL  MNTXEV 
PRINT  2  BLANKS 


•2  «  ANS? 


YES 


PRINT-FORMAT  (  FI )  -  (MC)MC) 
T)  X) EPI , F)  VPI,  F  )  SIG PI,  F) 

2H 


INC  (  F«  F+  3  ) 

NC  s  NC  I 
T:T»  OELT/2 
Xb  »  X 

X«  xe  SORT  (B/RHO) 
(OELT/2) 


SET 

(F.O)V.  2 

*  ANS 

NC  F 

-  MC  X:0  T*MC*DCLT 

SIBI 

*  SlOO 

El* 

(EPI,  F*  EP2,F 

)/2 

VI* 

(VPI,  Fl- VP2, 

F)/2 

SET 

(OAH*BIO.I) 

FN«A{0-EPl,(F-3) 


(l  +  B£TA»EPl,(F-S)) 


’  (l♦ALPHA«EPl.(F-3)} 

NC  >  NC  I 

SIOPI.  FtSlOPi,  (F-S) 

VPI,  F  «-  SIOPI.  F/C 
EPI,  Ft  SIO  PI,  F/e 


Uamilil 


NC*  N  I  »  ANS  5* 


TsT+  OELT/2 
Xb  *  X 

X  =  XB«-S0RT(^^)(5^)j  |3ET(3TlEXtB«  () 


PRINT-F0miAT{Fl)-(MC)NC)T)X)EPl,F)VPI,F)3I0l  ,F) 
2*1 


00  TO 
(.371  EX) 


OFPIs  0^^-|si6P2,Fs  0* 


GFPI  s  SI0P2,  F/EO  ♦ 

tlieXP(H  •  LM  (S1BP2,  F/eo)> 


SlOO  s  O  ^ 


NO 


6  ZERO«S]00/E0 


FCr*CPI.F(l«^MTA«EPl,n/(l-FALPHA  •  EPI.F) 
FEFP(*CP2,F<I*K7A»EP2^(I4-ALPMA»EPV^) 
FF  »  A  {02E  RO-FSF) 
FFPIsAIOFPI-FEFPl) 


FEl*  EI<l4-BETA»Ei)/(l4-ALPHA«  El) 
FI  .  A (0  ZERO- FED 


EMI  *  El 
VMI  *  VI 

EI*EPI.F<F0CLT(FF4FT)/(Z»B) 

VI  •  VPI,  (  F  ♦S)- (sjoo- SIOPI,  (FTS))/C 
-  (DCLTKFFPI«F7)/(4*C) 


SIO  P2,  F*0  ?  G  FPI  *  0 


NO 


GFPt*SIGP2,  F/EO  4- 

X*EXP(h*LOO  )) 


OF  s  0  1 - 1  SIOPI, Fs07 


OF*  SIOPI,  F/EO  ♦ 


k«exp(h«loo  ( 


OIMI  t  SIOI/EO  ♦ 
K*  EXP(k*LOO  ( 


SI  7 


I  *  ALPHA* EPI, F) 


(I«-BEJA*EI) 


(!«■  ALPHA*  El) 


FFPI 

A(SF»>i-FEFPr) 

FF 

A 

(OF-FEF) 

FIMI  * 

A 

(OIMI-FEJMI) 

«0  2«0  ?  — — ^  01*  0 


SIO  Ml 

*  SIO  2 

VM)  * 

V2 

EMI  t 

E2 

SIO  2 -SIB  Ml 


SIOMI 


39.2 


S0tll*O? 


CODE  2 


ro 

v-» 


9.  I 


CODE  3 


START 


f^:> 

ro 


BLOC  (EPl-EP7)/3 
(VPI-VP7)/3 
(SI0PI-ST0P7/3) 

FI  FORM  (9-IO>l-7)3*HO) 
LAST  (EPDVPDSIOPDl 


3  2 


BEGIN  ±  7 


|INC(F«  F»3)| 

_  i  vcj 

|siQPt.(F-3>=0?  I - •jo*  0  I 


{REAO(A)BEE]RHO)DELT)M)N) 

PRINT  LABELING  PRMT  BLANK 

PRIN T  FORMAT  (FI )(A)  BEE )RHO)OELT)MlN) 

PRINT  BLANK 

f  READ  ( H)  K  )EO )  BETA)  DELTA )  AL  PH  A ) 

I  PRINT  LABELING  PRINT  BLANK) 

\  PRINT  FORMAT(FIXH)K)EO)BrrA)OeLT)ALPHA) 
i^PRlNT  BLANK 
READ  ($tOO)BEeiN)  ANS) 

PRINT  labeling  print  BLANK 

PRINT-  FORMAT  (Ft)-(SKX))  BEGIN)  ANS)  B)  C) 

PRNT  2  BLANKS 

PRMT  LABELING  FOR  M  N  ETC 

PRINT  2  BLANKS 


_ IJ_ 

B iBEE • A 

SET  (Fs3)  C:S0RT(RH0*6) 

I  REAO'FORMAT  (FI)-(MC)  NC)  T)X) 
EPI.F)  VPI.F)3IQP1,F) 

I  SET  (3  71  EX  :  B  3  4  ) 

I  PRINT- FORMAT  (Fl)-(MC)NC)  T)  x: 
EPI.F)  VPl.F)  SIGPI,  F  ) 

I  z*i 


SET(3EX  :  871) 


INC(Fi 

F  +3) 

EPI,  F 

»0 

VPl,  F 

=  0 

SIO  PI , 

F:0 

r 


SET (FsO) 

REAO*  FORMAT  (Fl)-(MC)NC)T)X) 
EPDVPI)  SIGPI) 
WINT-FORMAT  (Ft)- (MC) NC)T)X) 
EPDVPl)  SIGPI) 

SIC  |:  SIGPI 

VI  s  VPl 

El  s  EPI 

1  «  I 

SET  (0AM  *  BIS  I) 


SI6P2.F7S»I 


|SET(3EXiBB  [  NC=N^ | 

NOl  4  2 

t  *  tNQ  f  Jii- 

[i  >  ANS  - [i  «E  »  I  \ 


PRINT  -  FORMAT  (FI) -{MC)NC)T)X) 
EPI.F)  VPI,F)SIGPI,  F) 

Z  :  I 


B7  I 


PRINT  2 
BLANKS 


0*S1GPI,  (F-3)/E0  ♦  K  • 
EXP(H«L00  (SIGPI, (F-3)/E0)) 


ITERATE  AGAIN 


|N0 

3  61 

SIGPI,  F-SIOMI 

<  OELTA 

SIGMI  1 

NO 

3.6 


3  32 


SRMlfO^ 


YES 
3  62 


SIO  PI.  Ft  SIO  PI.(F-3)-(0ELT/8)(A(G-€PI.F 


EPI,  FsSIGPI,  F/B 
VPl.  F. SIGPI,  F/C 


NCtNC^-l 
T«  T4-0ELT/2 
Xb  *X 


NC>  NC7I 

SIGPI,  FsSIGPI,  (F-3) 
VPl.F  f  SIGPI ,  F/C 
EPi.Fs  SIGPI,  F/B 


G*SIGPl,  F/EO  ♦ 
K*EXP(M«LOG 


MC>  M  ^ 


MC< 

PRINT  J 

MCf  1 

BLANKS 

MO 

B2I 

1 

MC 

:  M’ 

YES 

,  21  3 

PRINT 

BLANKS  1 

_ n 

1  Bldl 

T  t  T  f  BELT/ 2 

X0  *  X 

X  t  Xs  *s/b 

INC  (  FtFf  3 
NC«  NC4) 

/RNO  . 

1  ' 

SIGPI, FtO’ 


YES 


PRINT-FORMAT  (Fl)- 

(MC)  NC)T)X)EI)VI)SIGI) 
2  *  I 


T--  T+  OELT/2 
Xb*  X 

X.XB>SORT(-^)(fi^) 


_  3  93 

I  SET  (3  EX»BB  I  I 


— fn  NO  I  -  I  NO  t  - 1 

M?  ^  i-ANS’  2. 

YES _  I  1/^  <5'S'r 


I 


PRINT-  FORMAT  (Fl)-(MC)  NC)  T)  X )  EPI,  F )  SIGPI,  F  ) 

I  »  1 


YES 
I  133 


Z«ANS^I 


1  16  4 

.’'”1 

NO 

16 

SET  (GAM:  6  211 
EP2.  F  s  E2 

VP2,  F  :  V2 
S[0P2.  F  :  SIGI 

|nCsN^  I 

1 

•  1 

EP2,  F*  El 


Sts  P2,F.O  ?  Pl*0 


GFPfSlSPZ,  F/EO  * 
K-EXPtH«LOS  (- 


EO 

^IS 

1 

«F>St6Pt,  F/EO  4 

K*EXP(M«LOO(-^^i^^  9 


eiMI:  SJOl/EO 

K«  EXP(H»  LOO  >) 


FEFsEPl.F 

FEFP1«EP2,F 

FEIMIsEl 


(1*I>ETA»EPI.F) 
(If  alpha*  EPI.F) 

(  If  BCTA*EP2,F)  I 
(If  ALPHA*EP2,F)  \ 
(  I  f  BETA  •  El) 

(  If  ALPHA*EI} 


FFPIs  A  (OFPl-FEFPI) 
FF  *  A  (OF-  FEF) 
FIMI  ■  A  (0)MI  -FEIMI) 

I 


YES 


33  2 

*0 


33  I 

Jsioi'O? 

I 

|oi«si02/eof  I 

I  4TB  9  I 

K*EXP(MfLOO  ( 


(l^gCTA.EZ) 
(|fALPHA^E2) 
FI  »  A  (OI-FEI) 

I 


A 


A 


EPI,F 

El 

El 

E2 

VPl  ,F 

V  1 

VI 

V2 

SIGPI, F 

SIO  1 

SIGI 

$10  2 

A 


1 

NO 

37,1 

1  L  DELTA 

A  ^ 


(Sioif  S10P2^F>  C(VI-VPI,(F43)) 


^OELT^^iFlMI  f  FI)  ^  (FFPIf  Fl)^ 


SIG2-SieMI 

z.  OELTA 

^3 

1  1 

EMI  >0  F 

Slf  HI 

• 

E2i  EPI,  Ff 


2 

(SfOZ-SWPUF)  OCLr(FF»Fl) 


(SI02-S10I) 


at*B 

PELT  (FIMI  f-H) 

4*0 


"“i  s 

Isi02  »  0? 

V 


CODE  4 


INTERIOR  POINTS 


START 


M.l 


I  BLOC  (EP1-EP7)/3 
(VP1-VP7)/3 
( SIBP1-SIGP7/3 

I  FI  FORM(9-IO)I-7)SH>IO) 
LAST(EPI)  VPI)  S19PI) 


READ  (A)  BEE)  RHOI  OELT)M)N) 
PRNT  labeling  PRINT  BLANK 
PRKT  (A)  BEE)RH  )OELT)M)N) 
PRINT  BLANK 

'rEAO(H)K)E  )BETA)0ELTA)ALPHA) 
PRINT  LABELING  PRINT  BLANK 
PRINT  ( H )  K)  CO)  BETA)  DELTA )  ALPHA) 
PRINT  BLANK 
READ  (STGO)  BEGM)  ANS) 

PRifT  LABELN6  PRINT  BLANK 
PRINT-FORMAT  (Fi)-(S160)  BEGIN) 
ANS)B)C) 

PRINT  2  BLANKS 
PRINT  MN  LABELING 
PRINT  2  BLANKS 
B  s  BEE  •A  C:  V  RhO  /B 


SET  (  F.-O) 

READ-  FORMAT  (  Fi )  -  (MC)NC)T)  X )  EPI ) 
Vpl )  SIGpl  ) 

PRINT- FORMAT  ( F( )  •  {MC )  NOT)  X  )  E  PI) 
VPI  )  SIGPi) 

SIGI:  SIGPI 
VI  =  VPI 
El  :  E  PI 
Z  :  I 

SET  (GAM  :  B  19.1  ) 


READ-FORMAT  (  Fl )- MC)  NC ) T)  X  ) 
EPI,  F)  VPI,  F)  Stop,  F  ) 


ro 


MCt  1 

PRINT  2 

dC+l 

BLANKS 

NO 

B2I  1 

M  C  !  M  ? 

1 

_ 

2,  3 

PRIf^  2  BLANKS 

BI9  I 


T  *  T  +  OELT/2 

Xb-x 

X  «  Xb  Vb/rho  • 

INC  (F*  F-f  5)  * 

NC«  NC  ♦ 


PRINT-FORMAT(F|)-(MC)NC)T)X)EI)  VI)  SIGI  ) 
2  *1 


YES 

IS.3 


16  3 


I  SEr(GAM:B2.U)  | 
EP2.  F=  E2 
VP2,  Fs  V2 
SIGP2.  Fs  SI6  2 


2  :  ANS 

NO 

16.2 

NC: 

=  N  T 

PROGRAM  FOR  CODF  I  AND  INPUT  PARAMETERS  WITH 
OOUNDARY  DATA  FOR  A  SPECIAL  CASE 


FI 

START 


2.1 

3.1 


3.2 

3.3 


83. 4 

3.41 

3.42 

3.43 


3.5 


3.51 

3.52 


3.6 

3.61 

3.62 
3.71 

3.91 

3.92 


PROO  Tr-033  GROUiNu  SHOCK  STUDIES  CODE  1 
BL0C(EPI-EP7/3)VPl-VP7/3)5IGPI-SIGP7/3)% 

FORM  (<i-10)  1-7)5-1-10)5: 

LASTIEPDVPl  )SIGP1  )^ 

READ! A)BEE )RHO)DELT  )M)N)  ^ 

PRINT-FORMAT { F  1  )- 

CONT  <  A  BEE  RHO  DELT  > 

CONTC  M  N  >?  ENTER(PRINTB)* 

PRlNT-FOkMAT(Fl)-{ A)BEE)RHO)DELT)M)N)*ENTER(PRINTB)X 
READ! H)K)EO) BETA )OELTA) ALPHA) % 

PRINT-FORMAT ( F  1  )- 

CONT  <  H  K  EO  BETA  > 

CONT<DELTA  ALPHA  >*E NT E R { PR  I NTB ) * 

PRINT-FORMATIFl ) - ( H ) K ) EO ) BE T A ) DELT A ) ALPHA ) CENTER ( PR I N TB ) * 
READ(SIGO)BEGIN)ANS)^ 

PRINT-FORMAT { F  1  )- 

CONT  <  SIGO  BEGIN  ANS  B  C  >% 

ENTFRI PRINTB ) % 

IF (N<M)GOTO{ 2.2)X 
H  =  aEE*A5;  T  =  0:«X  =  05. 

Z  =  ANS:« 

SET{F=0) ^C=SORT{RHO*B  ) ?MC=0SNC=0* 

PR INT- FORMAT (FI )-( SIGO )BEG IN) ANS)B)C)* 

ENTER(PRINTR)?  ENTER! PRI NTB ) « 

PRINT-FORM.AT(Fl)- 

CONT  <  M  N  T  X  E  > 

CONT  <  V  SIGMA  >t 

ENTER(PRINTB)XENTFR(PRINTB)« 

IF(BEGIN)  IS  +  GOTO(4.  ID* 

SIGP1,F=SIG0*EP1 ,F=SIG0/B*VP1 ,F=-SIGO/C* 

2=1*  SET! 3. 71EX=B3.4) * 

PRINT-FORMAT! FI )-!MC)NC)T)X)EPI,F)VPl,F)SlGPI,F)% 

INC  !  F  =  F<-3)  * 

IF ! SIGPl  ,  ! F-3)  =0)G0T0(3.42)* 

G=SIGP1, !F-3) 

CONT  / E0+K*EXP ( H*L OG ! SIGPl,! F-3) /EO) ) *GOTO ! 3 . 43 ) * 

G  =  0* 

FN  =  A!G-EP1  ,  ( F-3 )  ! 1+BETA*EP1 , ! F-3)  )/(  1  +  ALPHA*EP1 , ( F-3)  )  ) * 

NC=NC+1%SIGPI ,F=SIGPi,!F-3)X 

VPl ,F  =  -S IGPl ,F /C*EP1 ,F  =  S IGPl  ,F/8% 

S1GM1=SI&P1 »F* 

IF!SIGP1,F  =O)G0TO!3.51)5: 

G=SI GPl ,F/EO+K*EXP ! H*LOG! S  IGP  1  ,F  /EO )  )*GOTO! 3. 52 ) * 

0  =  0* 

S1GP1,F  =  SIGP1 , ! F-3)-!0ELT/8) ! A(G-EP1,F( 1 +8 ET A*EP 1 , F  )  / 

CONT! 1+ALPHA*EP1  ,F)  )+FN)* 

EP1,F=SIGP1,F/B*VP1 ,F=-S IGPl ,F/C% 

IF!SIGM1=0)GOTO!3.62)* 

IF-ABS! ! SIGPl,F-SIGMl)/SIGMl<OELTA)GOTO!3.91 )*GOTO(3.5)1 
IF ! S IGPl ,F=0)G0T0! 3 .91 ) *GOTO ! 3.5)* 

PRINT-F0RMAT(F1)-!MC)NC)T)X)EP1,F)VP1,F)SIGP1,F)* 

Z=l*  GOTO!,3.71EX)* 

T=T+DELT/2*XB=X*X=XB+SQRT! B/RHO) *DELT/2* 

IF !NC=N)WI THIN! .005 )GOTO! 3.93)* 

21; 


1 

2 
3 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 
27 
2B 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


3.7  IF(Z=ANS)WITHIN( .005)0070(3.71)?  56 

3.64  Z=Z+1?  G0T0(B3.4)?  57 

3.93  SET( 3.71EX=e8. I )?GOTO( 3. 71 )?  58 

4.11  SET( 5eX=B7. L ) ?  59 

4.1  EP1,F=0?VP1 ,F=0?S1GP1,F=0?  60 

4.13  IFlNC^iNDWlTHlNC  .005)G0T0(4.  14)?  61 

4.2  IF(Z=ANS)WITH1N( .O05)G0TO(5.1)?  62 

4.3  Z=Z+1?  G0T0(B7.1)?  63 

4.14  SET( 5EX=B8. 1 ) ?  64 

5.1  PRlNT-FORyAKFn-lMONODXlEPl.FlVPltFlSlGPl.F)?  65 

Z=l?  G0TC(,5EX)9:  66 

B7. 1  1NC( F=F+3) ?NC=NC+1?T=T+DELT/2?XB=X?X=XB+SCRT (B/RHO) ♦DELT/2?  67 

G0T0(4,1)?  68 

B8.1  ENTER(PRINTB) ?ENTER(PR1NTB)?  69 

9.1  MC=1?  70 

11.1  SET(F=0)?  NC=RC?X=0?T=MC*OtLT?  71 

Z=ANS?  SET(GAM=B19. 1) ?  72 

SIGl  =  SIG03;El  =  (EPl,F+bP2,F)/2?Vl=(VPl,F+VP2tF)/2?  73 

11.2  1F(S1GP2,F  =0 )G0T0( 1 1 . 3) ?  74 

GFPL  =  SIGP2  ,F/E0<-K*EXP(  H*L0G(  S1GP2  ,F/E0)  )  ?GCT0(  11.35)?  75 

11.3  GFP1=0?  76 

11.35  1F(S1G0  =0)GOTO( 1 1.4)t  77 

GZER0=SlG0/6  04-K*EXP  (H*L0GI  SIGO/EO  )  )?G0T0(  11.5)?  78 

11.4  GZER0=0?  79 

11.5  FEF=EP1,F( ltBETA*EPl,F)/{l+ALPHA*EPl,F)?  80 

FEFP1=EP2,F(1+0ETA*EP2,F)/(1*ALPHA*EP2,F)?  81 

FF=A(GZER0-FEF)?  82 

FFP1=A(GFP1-FEFP1) ?  83 

12.1  FE l=t 1 { 1 +BETA*E1 ) / { 1 ♦ALPHA*E1 ) ?  84 

F1=A(GZER0-FE1 )?  85 

13.1  EMl  =  ei?Vf',l  =  Vl?  86 

E1=EP1,F+0ELT{FF+F1 )/(2*B)?  87 

V1=VP1 , {F+3)-(SlG0-SlGPl , ( F+3) ) /C-( DELT ) (FFP1+F1)/(4*C)?  88 

13.2  1F{EMI=0)G0T0( 13.3)?  89 

14.1  1F-ABS{ { El-EMl ) /EM1<DELTA)GOTO( 13.4 ) ?&0T0( 12.1 )?  90 

13.3  1F{£1=0)GCT0( 13.4) ?G0T0( 12. 1) ?  91 

13.4  1F{VM1=0)G0T0{ 13.5)?  92 

15.1  1F-ABS{ { V1-VR1)/VM1<OELTA)GOTO(15.2)?GOTO( 12.1)?  93 

13.5  IF{V1=0)G0T0{15.2)?G0T0{ 12.1)?  94 

15.2  EP1,F=E1?VP1,F=V1?  95 

16.2  IF(NC=N) G0T0( 16. 3) ?  96 

15.3  1F{ Z=ANS ) W1 THIN { .005)G0r0( 16- I ) ?  97 

15.4  Z=Z+1?  G0T0{B19.  1  )  S:  98 

16.1  PRINT-F0RKAT(F1)-{MC)NC)T)X)E1)V1)S1G1)?  99 

Z=l?  G0T0(C,GAM)?  100 

16.3  SET(GAM  =  821.  D?  G0T0{16.1)?  101 

B19.1  INC{F=F+3)?NC=NC+1?  102 

T=T+DELT/2?XB=X?X=XB+SQRT ( 8/RHO ) *DELT/2?G0T0 (30.1)?  103 

B21.1  IF(MC=M)WITHIN( .005)G0T0(21.3)?  104 

22.1  MC=MC  +  1?  ENTER(PRINTB)?  ENTER ( PR  I NT B ) ?  105 

22.2  IF(MC=-M)WITHIN(  .005)G0TO(22.3)?  GCTOdl.D?  106 

22.3  ANS=1?  GOTOdl.D?  107 

21.3  ENTER(PRINTB) ?  ENT ER ( PR INTB ) ?  GOTO(START)?  108 

30.1  SIG2=-(SIGP2,F  +  SIG1)/2?E2=(EP2,F+E1)Z2?  109 

V2= ( VP2tF*Vl ) /2?  110 

31.1  1F(SIGP2,F  =0)G0T0(31.2)?  Ill 

GFP1=SIGP2,F/E0+K*EXP(H*L0G( SIGP2,F/E0) ) ?G0T0 ( 31 . 3 ) ?  112 

31.2  GFP1=0?  113 


31.3  IF(SIGP1,F  =0  )GOTO(  31.A)  5i  114 

GF^SIGPl ,F/E0+K*EXP(H*L0G( SIGPl,F/e0) ) *G0T0 ( 3 1 . 5 ) *  115 

31.4  GF=0*  116 

31.5  IF(SIG1  =0)G0T0(31.6):«  117 

GIfU  =  SIGl/EO  +  K*EXP(  H*L0G(S IGl/eO) ) *GOTO ( 3 1 . 7 ) «  118 

31.6  GIM1=0«  119 

31.7  FEF=EP1,F  ( 1 +BE T A*E P I , F ) / ( 1 + ALPHA*EP 1 , F ) *  120 

FEFP1=EP2,F  ( U8ETA*EP2,F)/( 1+ALPHA#EP2,F)X  121 

FEIM1=E1  ( 1+BETA*E1 )/( 1+ALPHA*E1)%  122 

32.1  FFP1=A(GFP1-FEFP1) ?  F F = A ( GF-FE F ) ?  123 

F  IP1=A  (GIM-FEIMl )  %  124 

33.1  IF(SIG2  =0)G0T0(33.2 )«  125 

GI=SIG2/EO+K*EXP(H*LOG(SIG2/EO) ) XGGTO ( 33 . 3 ) *  126 

33.2  GI=0?  127 

33.3  FEI=E2( l+BETA*E2)/( 1+ALPHA*E2)X  128 

FI=A(GI-FEI ) ^  129 

34.1  SIG«1=SIG2XVMI=V2*EM1=E2%  130 

35.1  SIG2=(SIGl+SIGP2,F)/2-C(Vl-VPl,(F*3))/2  131 

C0NT-(DELT/4)((FIMl+FI)/2+(FFPl+FI)/2)*  132 

E2=EP1  ,F  +  (  SIG2-SIGP1,F )/B  +  OELT ( FF  +  F I ) / ( 2 *B ) *  133 

V2  =  V1+(SIG2-SIG1)/C  +  0ELT(FIM1+FI ) /( 4*0  *  134 

35.2  IF(SIGMl=0)G0T0(35.3)?  135 

36.1  IF-ABS(  (SIG2-SIGMn/SIGMl<0ELTA)G0T0(35.4)?G0T0(33.1)*  136 

35.3  IF ( S IG2=0 )GOTO ( 35. 4) XGOTO( 33. 1 )*  137 

35.4  IF(EM1=0)G0T0(35.5)X  138 

37.1  IF-ABS( (E2-EM1)/EM1<DELTA)G0T0(35.6)SG0T0( 33.1)S  139 

35.5  IF ( E2=0)G0T0( 35.6) ?GOTO( 33, I )?  140 

35.6  IF(VMI=0)G0r0(35,7)«  141 

38.1  IF-A0S{ (V2-VMl)/VMl<DELTA)G0T0(39,l)XG0T0t33.1)?  142 

35.7  IF(V2=0)G0TO(39.1)?G0T0(33,l)%  143 

39.1  EPI,F=E2?:E1=E22  VP  1 ,  F  =  V2tV  I  =  V2%  144 

SIGP1,F=SIG24SIGI=SIG2?G0T0( 16.2 )X  145 

2.2  PKINT<  N  IS  LESS  THAN  M  >*  146 

ENTERIPRINTBIX  ENTER ( PRI NTB ) «  GOTO(START)?  147 

END  GOTO{START)«  148 


CASE  1  10X10  DOMAIN  INPUT 

A  BEE  RHO 

90388944  07  23896900-02  15000000-03 
H  K  EO 

22200000  01  68600000  04  21400000  05 
SIGO  BEGIN  ANS 

50000000  02-50000000  00  10000000  01 

CASE  2  3X3  ZONE  INPUT 

A  BEE  RHO 

90388944  07  23896900-02  15000000-03 
H  K  EO 

22200000  01  68600000  04  21400000  05 
SIGO  BEGIN  ANS 

50000000  02-50000000  00  10000000  01 


PARAMETERS 

DELT  M  N 

20000000-03  10000000  02  10000000  02 
BETA  DELTA  ALPHA 

00000000  00  10000000-08  00000000  00 


PARAMETERS 

DELT  M  N 

20000000-03  30000000  01  30000000  01 
BETA  DELTA  ALPHA 

00000000  00  10000000-08  00000000  00 


26 


PROGRAM  FOR  CODE  2  AND  INPUT  PARAMETERS  WITH 
BOUNDARY  DATA  FOR  A  SPECIAL  CASE 


FI 

START 


3.1 


3.13 

3.14 

3.15 
9.  1 
11.1 


11.2 

11.3 

11.35 


1  1  A 


1 1  •  n 


11.5 


12.1 


13.1 


13.2 

14.1 

13.3 

13.4 

15.1 

13.5 


BLOC{EPl-EP7/3) VPl-VP7/3)SIGPl-SIGP7/3)« 

FORM (9- 10 ) 1-7) 5-1-101% 

LASKEPl  )VP1)SIGP1)% 

READ(A)BEE)RH0)0£LT)M)N)% 
print-format*  FI  )- 

CONT  <  A  BEE  RHO  DELT  > 

C0NT<  M  N  >%  ENTER(PRINTB)% 

PRINT-FORMAr(Fl 1 - ( A ) BEE ) RHO ) DE L T ) M) N ) %ENT ER ( PR INTB j % 

READ (H)K)EO) BETA) DELTA) ALPHA)? 

PRINT-FORMAT( FI  )- 

CONT  <  H  K  EO  BETA  > 

CONT<DELTA  ALPHA  >?ENTER( PRINTB) * 

PRINT-FORMAT! FI )-( H ) K) EOIBETA) DELT A ) ALPHA ) SEN TER ( PR  I NTB ) ? 
REAO(SIGO)BEGIN)ANS)? 

PRINT-FORMAKFl  )- 

CONT  <  SIGO  BEGIN  ANS  B  C  >% 

ENT£R(PRINTB)S 
8=BEE*A% 

SET(F=0)?C=SQRT(RH0*B  )? 

PRINT-FORMAT! FI )-( S IGO ) BEG  I N ) ANS ) B ) C ) S 
ENTERIPRINTB)?  ENTER! PRINTS )S 
PRINT-FORMAT!  FI  )- 

CONT  <  M  N  T  X  E  > 

CONT  <  V  SIGMA  >% 

ENTER! PRINTS) SENTE R I  PRINTS )%  M  =  N1 

READ-FORMAT ! F 1 ) - ! MC ) NC ) T ) X ) E P 1 , F ) VPl , F ) S I GP I , F ) ? 

IF!NC=N)WITHIN! .005)GOTOI9.l)? 

INC!F=F+3)?  G0T0I3.13)? 

MC=MC+1? 

SFT!F=0)%  NC=MC?X=OST=MC*OELT? 

2=ANS?  SET!GAM=B19.1)? 

SIG1=SIG0?E1=! EPl ,F+EP2,F) /2?V1=! VPl ,F+VP2t F )/2? 
iF!SIGP2,F  =0)G0T0! 11.3)? 

GFP1=SIGP2,F/E0+K*EXPIH*L0G!SIGP2,F/E0) )?G0T0! 11.35) ? 
GFP1=0? 

IF!SIG0  =0)G0T0! 11.4)? 

GZERO=SIGO/EO+K*EXP! H*LOGI SIGO/EO) )?G0T0! 11.5)? 

FEF=EPlfF! 1+BETA*EP1,F)/1 1+ALPHA*EP1,F)? 

FEFPl=EP2f F( 1  +  BETA*EP2.F )/ 1 1  +  ALPHA»EP2,  F ) ? 

FF=A!GZERO-FEF)? 

FFP1=A(GFP1-FEFP1) ? 

FEI=E1(1+BETA*E1)/! 1+ALPHA*E1)? 

FI=A(GZERO-FEI )? 

EM1=E1?VM1=V1? 

E1=EP1,F+DELT!FF+FI )/!2*8)? 

V1  =  VP1,  (F  +  3)-(SIG0-SIGPl,!F  +  3)  )/C-(  DELT)  (FFPH-FI)/(  4*0? 
IF(EH1=0)G0T0! 13.3)? 

IF-ABS! ! El-EMl ) /EM1<DELTA) GOTO ! 13.4)?G0T0! 12.1)? 
IF!E1=0)G0T0( 13.4) ?GOTO( 12.1)? 
IF(VM1=0)G0T0(13.5)?G0T0(I5.1)? 

IF-ABS ( ( V1-VM1)/VM1<DELTA)G0T0(15.2)?G0T0( 12.1)? 

IF! V1=0) GOTO (15.2) ?GOTOI 12.1)? 
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15.2  EPl  ,F  =  EnVPl  ,F  =  VU  56 

16.2  IF(MC=N)WITHIN{ .0O3)GOTO( 16.3)«  57 

15.3  IF{Z-ANS)WITHIN( .005)G0T0( 16.1)X  58 

15. A  Z=Z+1*  G0T0(6l9.1)t  59 

16.1  PRlNT-FOKKAr{Fi)-{MC)NC)T)X)ei)Vl)SIGl)*  60 

Z=l«  GOTOCOtGAM)?:  61 

16.3  SET{GAM  =  62l.l)^  G0T0(16.1)3:  62 

B19.1  INC(  F  =  F  +  3)?.NC=NC+1*  63 

T=T  +  0ELT/2:«X6=X:«X=XB  +  SQRT  (  B/ RHO  )  *OEL  T/2  %GOTO  (  30 . 1 )  *  64 

B21.1  IF(MC=M)WITHIN{ .005)G0T0(2l.3)«  MC=MC+1?  E NTER ( PR I NTB ) «  65 

ENTER{PRIN7B)t  I F { MC =M ) W I T H I N ( . 005 ) GOTO ( 22 . 3 ) %  G0T0<li.i)S  66 

21.3  ENTER(PRINTB)S:  EN  TER  (  PK 1 NTB )  *  GOTOISTART)?  67 

22.3  ANS=1%  G0T0(11.1)«  68 

30.1  SIG2={  S1GP2,F  +  S1G1) /2*E2=(EP2,F+E1)/2S>  69 

\/2={  VP2,F  +  V1  )/2t:  70 

31.1  IF(SIGP2,F  =0)G0T0( 31.2)%  71 

GFP1=S1GP2,F/E0+K*EXP( H*LOG( SIGP2,F/E0) ) XGOTO ( 31 . 3 ) %  72 

31.2  GFP1=0%  73 

31.3  IF(SIGP1,F  =0)G0T0(31.4)%  74 

GF  =  SIGP1  ,F/F0  +  K*EXP(H*L0G( S1GP1,F/E0) ) %GOTO ( 3 1 . 5 ) %  75 

31.4  GF=0%  76 

31.5  1F(S1G1  =0)G0T0{31-6)%  77 

GIM1=S1G1/E0+K*EXP{H*L0G(S1GI/E0) )XG0T0( 31.7)%  78 

31.6  GIM1=0%  79 

31.7  FEF=£P1,F  ( 1 +0E T A*£ P 1 , F ) / ( 1 + AL PHA*EP 1 t F ) %  80 

FEFP1  =  £P2,F  ( l+RETA*EP2,F)/(  l  +  ALPHA*£P2tF)%  81 

FEIM1=£1  ( l+6ETA*Fl )/( 1+ALPHA*E1)%  82 

32.1  FFP1=A(GFP1-FEFPI)%  FF=A( GF-FEF) %  83 

FUU=A(GIM1-FEIM1)  %  84 

33.1  IF(SIG2  =0)G0T0( 33.2 )%  85 

GI=SlG2/E0+K*exP(H*L0G(SlG2/E0) ) %GOTO ( 33 . 3 ) %  86 

33.2  GI=0%  87 

33.3  FEI=E2(l+BETA*E2)/{  1+ALPHA*E2)%  88 

FI=A{GI-FEI )%  89 

34.1  S1GM1=SIG2%VM1=V2%EM1=E2%  90 

35.1  SlG2=(SlGl  +  SlGP2,F)/2-C{Vl-VPl,(F  +  3)  )/2  91 

CONT-  { t)ELT/4  )  (  (  F  1  Ml+Fl  )  /  2  +  {  FFPl+Fl  )  /2)%  92 

E2=EPl,F+{ SIG2-SIGP1,F)/B+DELT{FF+FI )/( 2*B)%  93 

V2=Vl+{ SIG2-SIGI ) ZC+OELT { F I M 1 +F I ) / ( 4 *C ) %  94 

35.2  IF{ SIGM1=0 )GOTO ( 35. 3) %  95 

36.1  IF-ABS( 1 SIG2-SIGM1)/SIGM1<DELTA)G0T0C35.4)%G0T0(33.1)%  96 

35.3  IF{SIG2=0)G0T0{35.4)%G0T0{33.1)%  97 

35.4  IF ( EM1=0)G0T0( 35.5) %  98 

37.1  IF-ABS( ( E2-EM1 ) /EM1<DELTA)G0T0( 35.6) %GOTO( 33. 1 ) %  99 

35.5  IF{E2=0)G0T0(35.6)%G0T0(33.l)%  100 

35.6  IF (VM1=0)G0T0{35.7) %  101 

38.1  IF-ABS(  JV2-VMi)/V'Ml<0ELTA)G0T0(39.1)%G0T0(33.1)%  102 

35.7  IF(V2=0)G0T0(39.1)%G0T0{33.1)%  103 

39.1  EPl ,F=E2%E1  =  E2%  VP  1 , F  =  V2 %V I ^V2%  104 

S1GP1,F=SIG2%SIG1=S IG2%G0T0 ( 16. 2 ) %  105 

END  G0T0(START)%  106 
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H  K  EG  BETA  DELTA  ALPHA 

22200000  01  68600000  04  21400000  05  00000000  00  10000000-08  00000000  00 
SIGO  BEGIN  ANS 

50000000  02-50000000  00  10000000  01 

INPUT  60UN0ARY  DATA  ON  N=3  CHARACTERISTIC 

M  N  T  X  E  V  SIGMA 

300000  01  300000  01  bOOOOO-03  000000  00  451328-02-383640  02  500000  02 

300000  01  400000  01  700000-03  120000  01  412734-02-359208  02  477472  02 

300000  01  500000  01  800000-03  240001  01  379712-02-337303  02  456285  02 

300000  01  600000  01  900000-03  360001  01  351191-02-317622  02  436460  02 

300000  01  700000  01  100000-02  480002  01  326356-02-299895  02  417963  02 

CASE  2  INPUT  PARAMETERS 


A  BEE  RHO  DELT  M  N 

90388944  07  23896900-02  15000000-03  20000000-03  80000000  01  10000000  02 
H  K  EO  BETA  DELTA  ALPHA 

22200000  01  68600000  04  21400000  05  00000000  00  10000000-08  00000000  00 
SIGO  BEGIN  ANS 

5000C00C  02-50000000  00  10000000  01 

INPUT  BOUNDARY  DATA  ON  M=7  CHARACTERISTIC 

700000  01  700000  01  140000-02  000000  00  670435-02-459084  02  500000  02 

700000  01  800000  01  150000-02  120000  01  623943-02-438848  02  487296  02 

700000  01  900000  01  160000-02  240001  01  582185-02-419527  02  474410  02 

700000  01  ICOOOO  02  170000-02  360001  01  544555-02-401182  02  461515  02 
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PROGRAf'.  FOR  CODE  1  AND  INPUT  PARAMETERS  WITH 
BOUNDARY  DATA  FOR  A  SPECIAL  CASE 


PROB  TF-033  GROUND  SHOCK  STUDIES  CODE  3 

BLOC (EPl-EP7/3)VPl-VP7/3)SIGPl-SIGP7/3)% 

FI  FORMlO-10) 1-7)5-1-10)1 
LASTIEPl )VP1 )SIGP1)1 
START  REAO(A)REE)RHO)ueLTlM)N)l 

PRINT-FOKMATI FI )- 

CONT  <  A  BEE  RHO  UELT  > 

CONT<  M  U  >1  ENTERIPRiNTB)! 

PRINT-FORMAT (F 1 ) - ( A) BEE ) RHO ) DELT ) M ) N ) lENT ER ( P R  INTB ) 1 

REA0(H)K)F0)RETA)D£LTA)ALPHA)1 

PRINT-FORMATIFl )- 

CONT  <  H  K  EO  BETA  > 

CONTSDELTA  ALPmA  >1ENTER ( PRINTS) 1 

PRINT-FORMAT! FI ) - ( H ) K ) EO ) RE T A ) DEL T A ) ALPHA ) lEN TE R ( PR I NTB ) 1 
READ! SIGO) BEGIN ) ANS  )  1 
PRINT-F0RHAT(F1  )- 

CONT  <  SIGO  BEGIN  ANS  >1  ENTER ( PR  INTB ) * 

PRINT-FORMAT! FI )-<  S I GO ) B EG  I N ) ANS ) lENTER ! PR  I NTB ) lENTE R ( PR I N TB ) 1 
PRINT-FORMAT{F  I  )- 

CONT  <  M  N  T  X  E  > 

CONT  <  V  SIGMA  >1 

ENTEK(PRINT8)1ENTEK(PRINTS)1 

3.1  6=BEE«A1 
SeT(F=3)lC=S0RT(RH0*B  )% 

READ-FORMAT! FI ) - ( MC ) NC ) T ) X ) E P 1 , F ) VP 1 , F ) S  I GP 1 , F ) 1 
PRINT-FORMAT  1 F 1 ) - ( MC ) NC ) T ) X ) EP 1 , F ) VP  I , F ) S I GP 1 , F ) 1 
Z=ll  SET(3.71EX=B3.4)1 

3.2  IF!B£GIN)IS<-G0T0!4.  11)1 

B3.4  1NC!F=F+3)1 

IFISIGPl, !F-3)  =0)G0T0! 3.42)1 

3.41  G=SIGP1 , (F-3) 

CONT  /E0fK*EXPIH*L0G(SIGP1, (F-3)/E0) )1GOTO{3.43)1 

3.42  G=01 

3.43  FN=A!G-EP1,!F-3)!1+BETA*EP1,!F-3))/I1*ALPHA*EP1f(F-3)))1 
NC=NC+11SIGP1,F=SIGP1 , ! F-3 11 

VPl ,F=-SIGP1 ,F/C1EP1 ,F=SIGPl .F/Bl 

3.5  SIGM1  =  SIGP1  ,F1 

IFISIGPl, F  =0)G0T0!3.51)1 

G=SIGP1 ,F/E0+K*EXP! H*L0G! SlGPl ,F  /EO ) )1G0T0! 3.52)1 

3.51  G=01 

3.52  SIGPl .F^SIGPl , ! F-3 ) - ! DEL T/ B ) ( AIG-EPl ,F! 1 ♦BET A*EP 1 , F  )  / 
C0NT(1+ALPHA*EP1,F) )+FN)l 

EPl  ,F  =  SIGP1  ,F/81VPl ,F=-SIGP1 ,F/C1 

3.6  IF ! SIGM1=0)GOTo ! 3.62) 1 

3.61  IF-ABSI ! SIGPl ,F-SIGM1 ) /S I GM 1 <OEL T A ) GOTO ( 3. 9 1 )1G0T0! 3.5)1 

3.62  IF ( SIGPl ,F=0)G0T0( 3.91 ) IGOTOI 3.5) 1 

3.91  T=T+DELT/21XB=XlX=XB+SaRT! B/RH0)*DFLT/21 

3.92  IFINC=N)WITHINI . 005 ) GO TO ( 3 . 9 3 ) 1 

3.7  IF(Z=ANS)WITHIN( .005)G0T0I3.71)1 

3.64  Z=Z+11  G0T0(B3.4)1 

3.71  PRINT-FORMAT (FI )-!MC)NC)T)X)ePl,F)VPl,F)SIGPl,F)l 

Z=ll  GOTO!  ,3.71EX)  1 

3.93  SET! 3.71EX=B8. 1 ) IGOTO (3.71)1 

oci^DcA— or*x/-*> 
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4.1  INC(F  =  F  +  T)'i'  56 

EPl ,F=0^VPliF^0^SIGPl,F=0«  57 

4.13  IF(NC=N)WITHIN( .005)G0T0(4.14)X  58 

4.2  IF(Z=ANS)WITHIN( .n05)G0T0(5. 1)«  59 

4.3  Z=2+U  G0T0(R7.1)3!  60 

4.14  SET(btX=RR.i)%  61 

5.1  PRINT-F01<MAT(Fl)-(MC)NC)r)X)EPl,F)VPl,F)SIGPl,F)2  62 

Z  =  l?;  G0T0(,5EX):?  63 

87.1  NlC  =  i4C+l*T=T+oeLT/2%XB=X?X  =  XB  +  SQRT  ( B/RHO)  *DELT/2^  64 

G0T0(4.11)*  65 

B8.1  ENTER(PRIMTB)^ENTFR(PRINTB)«  GOTOdl.DX  66 

11. 1  SET!F=0)?  67 

READ-FORMAT  (FI  )-(  nC  )  iMC  )  r)X)EPl )  VPl  )  SIGPl  )*  68 

PRINT-F0KMAT(F1 )-( MC)NC )T)X)EP1 )VP1) SlGPl 69 
SIi;i  =  SIGPd  Vl  =  VPd  E1  =  EP1^  Z=l«  70 

SET  (GAM=B19.  1 )?;  71 

B19.1  Inc (F=F+3)*NC=NC+1«  72 

T=T+0ELT/2«XB=X%X=XB+SQRT( B/RHO) *DELT/2%  73 

30.1  SIb2=  (  SI(;P2  ,F  +  SIGl  ) /2%E2=(  EP2,F  +  E1  )/2^  74 

v2= ( vP2 t F + Vi ) /2?  75 

31.1  IF(SIGP2tF  =0)G0r0( 31.21 t  76 

GFP1=SIGP2,F/E0+K*EXP( H*LOG( SIGP2 ,F/EO) ) XGOTO ( 3 1 . 3 ) ?  77 

31.2  GFP1=0«  78 

31.3  IF(SIGP1,E  =0)G0T0( 31.4)%  79 

GF  =  SIGP1,F/E0  +  K*EXP(H*L0G(SIGP1,F/F0) ) XGOTO ( 3 1 . 5 ) «  80 

31.4  GF=0%  81 

31.5  1F(SIG1  =0)G0T0(31.6)%  82 

GIM1=SIG1/EO+K*EXP(H*LOG(SIG1/£0) )%G0T0( 31.7)%  83 

31.6  GIM1=0%  84 

31.7  FEF=EP1,E  ( 1 +fl E T A *E P 1 , F ) / ( 1 + ALP HA»EP 1 , F ) %  85 

FEFP1  =  £P2,F  {  l+BErA*EP2,F)/(  l-fALPHA*EP2,F)%  86 

F6IM1=£1  ( l+BETA*£l)/( 1+ALPHA*E1)%  87 

32.1  FFP1=A(GFP1-FEFP1)%  FF= A ( GF-FE F ) %  88 

FIM1=A{GIM1-FEIM1 )%  89 

33.1  IF(SIG2  =0)G0T0(33.2)%  90 

GI =SIG2/E0+K*EXP( H*L0G( SIG2/E0) ) %GOTO ( 33 . 3 ) *  91 

33.2  GI=0%  92 

33.3  FEI =E2 i 1 +BET A*E2 ) / ( 1 ♦ALPHA*E2 i *  93 

FI=A(GI-FEI )%  94 

34.1  Sir,MdSIG2%VKl  =  V2%EMl  =  E2%  95 

35.1  SIb2=(SIGl+SIGP2,F)/2-C{Vl-VPl,(F+3) )/2  96 

CONT- (  DELT/4  )  (  (  Flfll +F  1) /2  +  (  FFPl +F  1) /2  )4  97 

E2=EPlfF+( SIG2-SIGP1,F)/B+DELTIFF+FI )/(2*B)%  98 

V2=Vl  +  ( SIG2-SIG1 ) /C+OELT(F IMl+FI ) / (4*0  %  99 

35.2  IF(SIGM1=0)G0T0(35.3)%  100 

36.1  IF-ABS(  (  SIG2-SIGM1  )/S  I  GM  1<DELT  A  )  GO  TO  (  35 . 4  )  !«GOTO  (  33 .  1  )%  101 

35.4  IF( EM1=0)G0T0( 35.5) %  102 

37.1  IF-ABS( (E2-EMl)/EMl<0ELTA)GQT0I35.6)fG0T0(33. I)%  103 

35.5  IFIE2  =  0)G0T0(35.6).«G0T0(33.1)%  104 

35.6  1F(VM1=0)G0T0(35.7)%  105 

38.1  IF-ABS( ( V2-VM1)/VM1<UELTA)GOTOI39.1)%GOTO( 33.1)%  106 

35.7  IF(V2=0)G0T0«39.1)%G0T0{33.1)%  107 

35.3  IF ( SIG2=0)G0T0( 35.4 ) %G0T0( 33. 1) %  108 

39.1  EPltF=El%  E1=E2%  VP1.F=V1*  V1=V2*  109 

SIGP1,F=SIG1%  SIG1=SIG2%  G0T0(16.2)%  110 

B21.1  IF { HC=H i Wi THIN! . 005 )GOTO 1 2 i . 3 ) %  111 

22.1  MC=MCF1%  ENTER(PRINTB)%  ENTER! PRINTS)*  112 

22.2  IF(MC=M)W1THIN(  .005)GOTO(22.3)*  GOTOdl.D*  113 
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22.3  ^NS=lt  GOTO(ll.l)'^  114 

21.3  ENrER( PRINTB ) ^  ENTER( PR INTB ) ?  GOTOJSTART)*  115 

16.2  IF(NC  =  iN)GOTO(  16.3)  116 

15.3  I F  (  Z=AiNS)WI  THIN(  .005)  GOTO)  16. 1)%  117 

15.4  Z  =  Z+1*  G0r0(B10.  1) 118 

16.1  PHI'MT-FnRMAT(Fl)-(MCINC)T)X)£l)Vl)SIGl)X  119 

Z=U  G0T0(0,GAM)?  120 

16.3  SE r( GAM=B21 . 1 ) *  EP2,F=E2*  VP2tF=V2«  S I GP2 , F= S I G2^  121 

G0T0(16.1)%  122 

END  GOTO  (  START  )  123 


INPUT  PARAMETERS 


A  BEE  RHO  DELT  M 

90388944  07  23896900-02  15000000-03  20000000-03  30000000  01 
H  K  EO  BETA  DELTA 

22200000  01  68600000  04  21400000  05  00000000  00  10000000-08 
SIGO  BEGIN  ANS 

50000000  02-50000000  00  10000000  01 

INPUT  BOUNDARY  DATA  ON  N=3  CHARACTERISTIC 

000000  00  300000  01  300000-03  360001  01  183368-02-220043  02  396079  02 

100000  01  300000  01  400U00-03  240001  01  260629-02-275758  02  438797  02 

200000  01  300000  01  500000-03  120000  01  350227-02-330727  02  473202  02 

300000  01  3U0000  01  600000-03  000000  00  451328-02-383640  02  500000  02 


N 

10000000  02 
ALPHA 

00000000  00 
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FI 

START 


1.1 

11.1 


B19.  1 

16.2 
16.  J 
15.4 
16.  1 

16.3 
B21 . 1 
22.1 
22.2 

22.3 
2  1.3 
30.  1 

31.1 

31.2 

31.3 

31.4 

31.5 


PROGRAM  FOR  CODE  4  AND  INPUT  PARAMETERS  WITH 
BOUNDARY  DATA  FOR  A  SPECIAL  CASE 


PROB  TF-033  GROUND  SHOCK  STUDIES  CODE  4 

BLOC  itPl-EPT/'i)  VPI-VP7/3)  SlGPl-SlGP7/3)% 

F0RH(9-10)  l-7)5-l-iO)S; 

LAST (EPl ) VPl  )SIGPI )  % 

READ(  A)BlE  )RHO)DELT  )M)N);g 
PRINT-FOKf,AT(  FI  )- 

CONT  <  A  BEE  RHO  DELT  > 

CONT<  M  N  >%  ENTER(PRINTB)? 

PRINT-FORMAT  (  FI  )-(  A  )REE  )  RHO  )DELT)M)N)35ENTFR(PRINTP):S 
READ (H)K)EO) BETA) DELT A) ALPHA)* 

PRINT-FORMAT ( FI  )- 

CCNT  <  H  K  EO  BETA  > 

CONT<DELTA  ALPHA  >*ENTER ( PR I NTB ) * 

PRINT-FORMAT ( F 1 ) - { H ) K ) EO ) RET A ) DE LT A ) ALPHA) CENTER (PRINTB)* 
READ! SI  GO) BEG  IN) AMS  )* 

PRINT-FORMAT! FI  )- 

CONT  <  SICO  BEGIN  ANS  B  C  >* 

ENTEIMPRINTH)* 

R=BEE*A*  C=SURT( RHO*B)* 

PRINT-FORMAT! FI )-( S I  GO ) BEG  I N ) ANS  )  B ) C  )  * 

ENTER!PR1NTB)*  ENTER ! PR INTB ) * 

PRINT-rORMAT!Fl)- 

COl^/T  <  M  N  T  X  E  > 

CONT  <  V  SIGMA  >* 

ENTERIPRIMTH) *LNTER!PRINT8)* 

SET! F=3)* 

READ-FORMAT ! F 1)-! MC )NC)T )X )EP1,F )VP1,F)SIGP1,F)* 

IF!NC  =  N)WITHIN!  .005)G0T0!11.1)*  INC!F=F  +  3)*  GOTOll.D* 

SE  I  !F  =  0)i;  RFAU-FOkMAT!Fl)-(MC)NC)T)X)EPl  )VP1)SIGP1  )* 
PRINT-FORMAT! FI ) -  I MC ) NC ) T ) X ) EP 1 ) VP  I ) S I GP 1 ) ?  SE T I G AM  =  B 1 9 . I ) * 
SIG1=SIGP1*  V1=VP1*  EI  =  EPU  Z=l*  G0T0!B19.1)* 

INC  !F  =  F  +  3)?:NC  =  NC  +  1* 

T= l+0ELT/2*X0  =  X*X  =  XB  +  S0RT( B/RHO) *DELT/2*G0T0! 30. 1  )* 
IF!NC=N)WITHIN!.005)GCT0(16.3)* 

IF!Z  =  ANS)WITHIN! . 00 5 ) GO TO !  I  6 . 1 ) * 

Z=Z+1*  G0T0!B19.1)* 

PRINT-FORMAT ! FI )-! MC INC )T)X)EI)V1)SIG1)* 

1=1%  G0T0!0,GAM)* 

SETIGAM=B21.1)*  EP2,F=E2*  VP2,F=V2*  S ICP2 , F = S I G2*  G0T0!16.1)* 
IF!MC  =  M)WITHIN!  . 005 ) GO  TO ! 2 1 . 3 ) * 

MC=MC  +  1*  ENTERIPRINTB)*  ENTER ( PR  I NT B ) * 

IE  !MC  =  M) wl THIN! .005)G0T0 !22. 3) *  GOTO! 11. 1 ) * 

ANS=1*  GOTOIll.D* 

ENTER!PRINTB)*  ENTE R ! PR INTB ) *  GOTO!START)* 

SIG2=!SIGP2,F+SIG1 ) /2*E2 = ! E P2 , F +E 1 ) /2 * 

V2=! VP2,F+V1 )/2* 

IF!SIGP2,F  =0)G0T0! 31.2)* 

GEP1=SIGP2,E/E0+K*EXP(H*L0G! SIGP2.F/E0) ) *GOTO ! 3 1 . 3 ) * 

GFP1=0* 

1F!SIGP1,E  =0)G0T0!31.4)* 

GF=SIGP1,E/E0+K*EXPIH*L0G( SIGP1,F/E0) ) *GOTO ! 3 1 . 5 ) * 

GF  =  0* 

IF!SIG1  =0)G0T0! 31.6)* 

GIM1  =  SIG1/E0  +  K*EXP! H*L0G(SIG1/E0) ) *GOTO ! 31 . 7 ) * 
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31.6  GIf11=0?  56 

31.7  FGF  =  CP1,F  (  1 +tiFTA*EPl  ,F)  /  (  1+ALPHA*EP1  ,F )  *  57 

FEFP1  =  EP2,F  (  l+BETA*EP2,F)/( 1+ALPHA*EP2,F)X  58 

FEIM1=E1  ( 1+BETA*E1 ) /( 1+ALPHA*E1 )%  59 

32.1  FFP1=A(GFP1-FEFP1):«  F  F  =  A  ( GF-FE  F )  .«  60 

FlMl=A(GIMl-FEIKl)it  61 

33.1  IF(SIG2  =0)G0T0( 33.2)^  62 

G1=SIG2/EO+K*EXP(M*LOG(SIG2/EO) ) ?G0T0 1 33 . 3 ) «  63 

33.2  GI=0^  64 

33.3  FEl  =E2  (  1+BETA*E2 )  /  (  1 +ALPHA*E2  li’  65 

FI=A(GI-FEI 66 

34.1  SIGMl=SIG2%VMl=V2iEMl=E2^  67 

35.1  S1G2=( SIG1+SIGP2,F ) /2-C ( Vl-VPl, (F+3) ) /2  68 

C0iNlT-(L)ELT/4)((FIMl+FI)/2+(FFPl+Fl)/2)%  69 

E2=EP1 ,  F  +  (  S  IG2-S  IGP  1,F  )  /  B  +  OELT  {  FF  +  FI  )  /  (  2*B  )  9!  70 

V2=V1+ ( S IG2-SIG1 ) /C+UELT( F IMl+F 1 ) / (4*C)«  71 

35.2  IF  (  SIGM1  =  0)  G0TU(  35. 3)  9;G0T0(  36. 1  ) «  72 

35.3  IF ( SIG2=0) G0T0( 35.4 ) ^GOTO ( 33. 1 ) *  73 

35.4  IF(  EMl  =  0)G0r0(  35.5)  9:G0T0(  37.  1)%  74 

35.5  IF ( E2=0) GOTO( 35.6 ) XGOTO( 33. 1 ) ^  75 

36.1  1F-ABS( ( 5IG2-S1GM1 ) /SIGK1<0ELTA )GOTO( 35.4) ?GOTO( 33. 1 ) «  76 

37.1  IF-ABSl (E2-EM1)/EM1<DELTA)G0T0(35.6)XG0T0(33.1)«  77 

35.6  IF ( VM1=0 )G0T0( 35. 7) «G0T0 ( 38. 1 ) *  78 

35.7  IF  ( V2  =  0 )  G0T0(  39. 1  )  :«G0T01  33.  1)  9!  79 

38.1  IF-A8S(  ( V2-VrU  ) /VM1<DELTA) GOTO! 39- 1 ) SGOTOl 33. 1 ) «  80 

39.1  EP1,F  =  EUE1  =  E2'«  VP  1  ,  F  =  V 1  i  V 1  =  v2  91  81 

S1GP1,F=SIG1%SIG1=SIG2^G0T0(16.2)?  82 

END  GOTO  (  START  )  9:  83 


INPUT  PARAMETERS 

A  BEE  RHO  DELT  M  N 

90388944  07  23896900-02  15000000-03  20000000-03  70000000  01  10000000  02 
H  K  EO  BETA  DELTA  ALPHA 

22200000  01  68600000  04  21^00000  05  00000000  00  10000000-08  00000000  00 
SIGO  8EGIN  ANS 

50000000  02-50000000  00  10000000  01 

INPUT  BOUNDARY  DATA  ON  M=3  CHARACTERISTIC 

300000  01  700000  01  100000-02  480002  01  326356-02-299895  02  417963  02 

300000  01  800000  01  110000-02  600002  01  304572-02-2 8 3 8 8 2  02  400730  02 

300000  01  900000  01  120000-02  720003  01  285339-02-269371  02  384683  02 

300000  01  100000  02  130000-02  840003  01  268255-02-256182  02  369738  02 

INPUT  80UN0ARY  DATA  ON  N=7  CHARACTERISTIC 

400000  01  700000  01  110000-02  360001  01  402989-02-342221  02  444541  02 

500000  01  700000  01  120000-02  240001  01  486855-02-383498  02  467066  02 

600000  01  700000  01  130000-02  120000  01  576599-02-422768  02  485524  02 

700000  01  700000  01  140000-02  000000  00  670435-02-459084  02  500000  02 


34 


V.  RESULTS 


An  example  is  computed  on  a  grid  of  maximum  m  =  maximum  n  =  10^  both  by  use 
of  code  1  and  by  use  of  the  set  of  codes  1,2,3,^  (Figure  6). 


Code  1  uses  data  along  the  wave  front  and  boundary;  code  2  uses  data  along 
the  smallest  m  of  its  zone  and  along  the  boundary;  code  5  uses  wave  front  data  and 
data  along  the  smallest  n  of  its  zone;  code  U  uses  data  along  the  smallest  m  and 
smallest  n  of  its  zone. 
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CCMPLETE  10X10  DOMAIN 
USING  CODE  1. 
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COMPLETE  10X10  DCMAIN 
USING  CODES  1,  2,  5,  4. 
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